KEK-CP-126 
June 2002 



Matching between matrix elements and parton showers using a 
Leading-Log subtraction method in NLO-QCD 



Y. Kurihara 

. . . High Energy Accelerator Research Organization, 

, Oho 1-1, Tsukuba, Ibaraki 305-0801, Japan 

o 

(N 

Abstract 

A new method to construct event-generators based on next-to-leading order QCD matrix- 
' elements and leading-logarithmic parton showers is proposed. Matrix elements of loop di- 

agrams as well as those of a tree level can be generated using an automatic system. A 
soft/coUinear singularity is treated using a leading- log subtraction method. Higher order 
re-summation of the soft/collinear correction by the parton shower method is combined 
with the NLO matrix-element without any double-counting in this method. 
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' 1 Introduction 

o ■ 

O 



The standard model(SM) of electromagnetic and strong interactions has been well established 
through precision measurements of a large variety of observables in high-energy experiments 
However, there is still one missing part of the SM to be confirmed by experiments, i.e. the Higgs 
boson. In order to search for those missing parts and also for direct signals beyond the SM, the 



Q I LHC experiments [0 are under construction at CERN. Though LHC has employed a proton- 
• proton colliding machine due to a view of the discovery potential for new (heavy) particles, 
a capability for precision measurements must also be seriously considered. For extracting 
\ physically meaningful information from a large amount of experimental data, including huge 
■ QCD backgrounds, a deep understanding of the behavior of QCD backgrounds is essentially 
^ important. 

The prediction power of a lowest calculation of QCD processes is very limited due to a large 
coupling constant of the strong interaction and an ambiguity of the energy scale. Predictions 
based on the NLO (or more) matrix elements combined with some all-order summation is 
desirable. Much effort to calculate NNLO corrections!^ is being made for LHC experiments. 
Moreover, in order to match with the experimental requirements for precision measurements, a 
simple correction factor (so-called iC-factor) for the lowest calculation is clearly not sufficient. 
One must construct event generators including NLO matrix elements with higher order re- 
summation even for processes with multi-particle final states. When one tries to do this task, 
one may encounter the following difficulties: 

• the number of diagrams contributing to processes with many particles in a final state is 
very large, 

• any numerical instability in the treatment of a collinear singularity must be taken care 
of, 

• a careful treatment of matching between matrix elements and the re-summation part 
must be required to avoid double counting, and 



• any negative weight of events in some part of phase space is not avoidable. 



In order to overcome these difficulties, lots of work has already been done|^, Q . In this report, we 
propose a new method to construct event-generators with next-to-leading order QCD matrix- 
elements and leading-logarithmic parton showers. Our solution for the above problems, except 
for the last one, will be given in this report. A treatment of the negative weight will be presented 
elsewhere. 



2 NLO matrix elements 

A calculation of the matrix elements of hard-scattering processes is not such a simple task 
when multiple partons are produced. At the tree level, the GRACE systemj^], which is an 
automatic system for generating Feynman diagrams and a FORTRAN source-code to evaluate the 
amplitude, is established. The system can in principle treat any number of external particles, 
and has been used for up to six fermionsQ in the final state within a practical CPU time. In the 
GRACE system, matrix elements are calculated numerically using a CHANEL library based on a 
helicity amplitude. CHANEL contains routines to evaluate such things as: wave-functions/spinors 
at external states, interaction vertices, and particle propagator. 

For the loop amplitude, effective vertices of two- and three-point functions have been 
prepared^ and implemented in the system. Here, we employed the MS schemeQ for renor- 



malization and dimensional regularization|10] for treating soft/collinear singularities. For more 
than three-point diagrams, an extended CHANEL is prepared. It can evaluate a fermion current 
including any number of gamma matrices. An amplitude of loop diagrams will be numerated 
based on the fermion current connected by boson propagators. When a fermion loop is included 
in the diagram, numeration will be performed after taking its trace. A bosonic loop can also 
be treated by the CHANEL. 

For a loop integral of four-point functions, general formulae of the tensor integral are pre- 



pared and implemented in a library. For more than four-point functions, a reduction formula] 11] 
is used in the system. 



3 Event generation with a parton shower 
3.1 Conventional method 

The calculations of cross sections for proton-(anti) proton collisions are usually performed as 
follows m, 11: 

1. matrix elements of a given process at a parton level are prepared, 

2. a probability to observe a parton in a proton at given energy scale and momentum fraction 
is obtained from some parton distribution function(PDF) , 

3. a numerical integration is performed in phase-space with some initial energy distribution 
in the PDF, 

4. for initial-state partons, a (backward) parton shower [1l5|] is applied from the given (high) 
energy scale to a low energy-scale of partons, 

5. for final-state partons, a parton shower is also applied from a given energy scale to close 
to the hadron energy-scale, 

6. some non-perturbative effects according to string fragmentation or color clustering are 
taken into account, and 

7. physical hadrons are formed based on some hadronization model. 



In this method, the parton shower(PS) is used in a supplementary way, just for generating 
multi-partons with a finite transverse momentum. 



3.2 An x-deterministic PS 



The parton shower is a method to solve a DGLAP evolution equation 1 16 1 using a Monte Carlo 
method. The PDF's are also obtained by solving the DGLAP equation with experimental 
data-fitting. When an initial distribution of partons is given at some energy scale, the PS can 
reproduce a consistent result as the PDF. 

The PDF can give the weight of partons when the momentum fraction (x) and the energy 
scale are given by the users. On the other hand, in the forward evolution scheme of the PS, 
the momentum fraction of a parton (x) is determined only after evolution takes place. This 
method is very inefficient, for instance, for a narrow-resonance production. In order to cure this 
inefficiency, we have developed an x-deterministic PS. In a Monte-Carlo procedure in the PS, the 

evolution is controlled by the Sudakov form-factor, which gives a non-branching probability 
of partons, and the x determination is done independently according to a splitting function. It 
is not necessary to determine the x value at each branching. After preforming all branching 
procedures using the Sudakov form-factor, the x value at each branching is determined to give 
a total X, being a given value from outside of the PS. In this PS, each event has a different 
weight according to the splitting functions and initial distribution of the PDF. However, if we 
employ an appropriate transformation of the input random numbers, we can expect a sufficient 
efficiency with a weight close to unity. 

For singlet partons, we employ an evolution scheme based on momentum distributions 
rather than those on particle-number distributions. This scheme has been proposed by Tanaka 
and Munehisa|17] and shows numerically stable results with high efficiency. 

A numerical test of the x-deterministic PS is done by comparing results with the PDF, 
Cteq5L[Q. The initial distributions of partons are taken from Cteq5L at the energy scale of 
the 5-quark mass. Then, after evolution to an energy scale of 100 GeV, the x distributions from 
the PS are compared to those in Cteq5L. One can see in Fig.|| the PS reproduced x distributions 
in Cteq5L for both of singlet and non-singlet partons. 

Our procedure to generate QCD events involves the following three steps, : 

2'. numerical integration is performed in the phase-space including x integration, 

3'. the probability for the existence of each parton in the proton at some low-energy scale is 
obtained from the PDF, and 

4'. the x-deterministic PS performs evolution from a low energy scale to the energy scale of 
a hard scattering while generating multi-partons with a finite transverse momentum, 

instead of steps 2, 3, and 4 in the conventional procedure. 



4 Real emission 

4.1 Soft/Collinear treatment 

At first, the conventional method of phase-space slicing is used to treat soft and collinear 
singularities. It is explained for the case of an initial-state radiation in this report. The final- 
state radiation can be treated in a similar way. Let's consider two colored partons, whose 
momenta are and scattered into the A^-body final state of colorless particles as a Born 
process. For the radiative correction of this process, we must treat processes with one additional 
colored parton, whose momentum is ^ emitted in addition to the Born process. The matrix 
elements of the real emission processes must be integrated in (A^-l-l)-body phase-space, ^tv+I' 



in a d-dimensional space-time. The space-time dimension is set to be d = 4:+2eiji in this report. 
A coUinear region in the {N + l)-body phase-space of the final particles in the d-dimensional 
space-time is defined as 
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where is some cut-off value. It must be sufficiently small so as not to be observed exper- 
imentally. A final result must be independent of this value. In the collinear region, matrix 
elements can be approximated as 
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where p(^)(a;) is a splitting function for a given parton splitting at a leading-logarithmic order, 
and k± is the transverse momentum of the radiated parton and fi the energy scale of the 
splitting. The CM-energy of the Born process is s = = sx, where s = (pi +^2)^- Here, we 
can set d to be four in the matrix element of the Born process, since it has no IR singularity. 
In the same approximation, the phase space is expressed as 
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The total cross section of real-radiation processes in the collinear region can be obtained as 
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where ^vis means a visible region of the phase space, such as 

^vis ~ ^full ^coll- 

Then, the collinear cross section can be obtained after integration as 
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where 
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is the Born cross-section at a CM energy of 



xs 



in four-dimcnsional space-time, and 



L = In [s/ ji^), is a factorization energy scale, and /c is a color factor of the given branching. 
The IR-divergent terms in the coUinear cross section can be canceled out after summing up 
terms from virtual corrections. Even after combining with virtual corrections, there still exists 
the IR-divergent term, (pix^em) = ^7^^^^"^^^)- This term is thrown away by hand, since it is 
counted in the PDF or PS on the initial partons. The collinear remnant in the visible region 
{(^vis) gives large negative values when the threshold energy {Ql) is sufficiently small. This 
large negative number will be compensated by the large positive cross section of real-emission 
processes. We propose to combine cr^j^ with exact matrix elements of the real emission processes 
before the phase-scape integration. It will be discussed later in this report. 



4.2 Visible jet cross-section 

The last part of the NLO cross-section calculation is a real emission of the additional parton 
into the visible region with exact matrix elements. Those matrix elements of (A^+l)-parton 
production at the final state in four-dimensional space-time can be automatically generated 
using the GRACE system. Cross sections are obtained by integrating those matrix elements 
under the four-dimensional phase-space as 
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The space-time dimension is set to be four {em = 0) hereafter, since there is no IR-divergence 
in Qvis- When the threshold energy (Q^) is set to be sufficiently small, this cross section can 
be larger than the Born cross section due to the largeness of the coupling constant A 
parton-level calculation has no problem, except for this largeness of the higher order correction. 
However, if one tries to calculate the cross sections of a proton (anti-)proton collision, one has 
to combine the matrix elements with the PDF or PS to construct a proton from partons. The 
PDF and PS include leading-log(LL) terms of the initial-state parton emission. If one combines 
the matrix element with the PDF or PS very naively, one cannot avoid double-counting of these 
LL-terms. Our proposal to solve this problem is to subtract the LL-terms from the exact matrix 
elements as 
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The second term of the integrand is the LL-approximation of the real-emission matrix-elements 
under the collinear condition. There is nothing but V^js' in the last subsection. This 'LL sub- 
traction method' can cure not only the double-counting of LL-terms in the real-emission part, 
but also any large cancellation between the collinear remnant and the real-emission correction. 
When the threshold value is set to be sufficiently small, the integrand in aLLsub is very close 
to zero, because the LL-approximation is very precise around the collinear region. Then, the 
result of integration in a^Lsub is independent of the threshold value of Q^. 

This subtraction may distort the experimental observables of additional jet distributions. 
However, a subtracted LL-part will be recovered after adding the PS applied in the Born 
process. The exact matrix element of real-radiation processes gives only non-logarithmic terms 



to the visible distributions in the LL-subtraction method. The relation between the PS and 
the LL-subtraction method is as follows. The LL-approximation cross section in (JLLsub is 
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where Q^ax is the maximum value of the transverse momentum. V^^^x) is just the first term 
of the all-order re-summation of the leading logarithmic terms in the PS. Actually, the Sudakov 
form factor can be obtained from T)^^\x) as: 



n(QL.,Qc) = exp^-£vWix)dxy 



However, there is one essential difference: the upper bound of k'^ integration. In the PS, this 
integration is performed up to the energy scale determined by the Born process. Then, if aLL 
is subtracted in the full phase space of the (A^-l-l)-particle final state, some part of the LL cross 
section cannot be recovered from the Born cross-section with the PS. Then, an appropriate 
restriction on the phase-space integration must be applied, such as 
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where Q\ is the energy scale of the Born process. 



4.3 A test of the LL-subtraction method 

The LL-subtraction method is numerically tested for a process of uu /^^^^ in proton-anti 
proton collision at the CM energy of 2 TeV. Here only non-singlet u-quark is used in this test 
this test to avoid an additional complexity. Matrix elements of the Born process and the real 
radiation process {uu — > fi^^i~gluon) are generated by GRACE. Numerical integration is done 
using a BASES |jl^ system. A non-singlet w-quark distribution at the energy scale of 4.6 GeV is 
taken from CTEQ5L. The parton evolution from this energy scale to that of hard-scattering, i.e. 
an invariant mass of the muon-pair (s^^), is done using the x-deterministic PS for both of the 
Born and radiative processes. An additional cut of ^/s^ > 40 GeV is also applied. 

The distributions of transverse momenta of gluons are shown in Figure 2. For the Born 
process, the largest kj_ from the PS is filled to histograms if it has a k± greater than 1 GeV. 
For the radiative process, that from matrix elements is filled when it has passed the same cut 
as mentioned above. If one does not take care of double-counting of the LL-terms in the matrix 
elements and the PS, the distributions of da exact / dk_\_ show larger values compared with the 
PS in the Born process, as shown in the left histogram of Figure 2. In order to avoid double- 
counting, we required that k±_ of the matrix elements be greater than those of any gluons from 
the PS in the radiative process. This double-counting rejection has a good agreement of the k_\_ 
distribution around the low-fej_ region, as shown in the middle histogram. However, the PS can 
still emit an insufficient number of gluons around a very high-fe^ region. If the upper bound of 
the k\ integration is set to be s^^, the distribution from exact matrix elements shows a good 
agreement with that from the PS, as shown in right histogram. 

The distributions from aLLsub combined with those of the PS on the Born process are 
compared with those from exact matrix-elements with double-counting rejection. In both of 
the gluon energy and the transverse-momentum distributions, the LL-subtraction method can 
give consistent results with those of the exact matrix elements, as shown in Figure 3. 



5 Conclusions 



A new method to construct event-generators based on next-to-leading order QCD matrix- 
elements and an x-deterministic parton shower is proposed. Matrix elements of loop diagrams 
as well as those of a tree level can be generated by the GRACE system. A soft /collinear singularity 
is treated using the leading-log subtraction method. It has been demonstrated that the LL- 
subtraction method can give good agreement with the exact matrix elements without any 
double-counting problem. 
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Figure 1: Momentum fraction distributions from an x-deterministic parton shower and Cteq5L 
parton distribution functions. The distributions of u-quarks (upper left), iZ-quarks (upper 
right), and gluons (lower left) are shown. 




50 100 

P't (GeV) 




50 100 

p'r (GeV) 




50 100 

p't (GeV) 



Figure 2: Transverse momentum distribution of gluons. The distributions from tlie PS applied 
to the Born process are shown by * commonly in three histograms. Those distributions from 
o'exact (left), aexact with double-count rejection (middle), and (Texact with double-count rejection 
and the k\ restriction are compared with the PS. 
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Figure 3: Energy and transverse momentum distributions of gluons. The * shows the dis- 
tributions from the PS on the Born process, sohd histograms from Gexact with double-count 
rejection, and circles from aiLsub combined with those from the PS on the Born process. 



